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. Abstract. ThermaUzation of quark-gluon matter is studied with a transport 

equation that includes contributions of 2-to-2 and 3-to-3 elastic scatterings. 
Thermalization time is related to the squared amplitudes for the elastic scatterings 
that are calculated in perturbative QCD. 
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For a nonequilibrium partonic system, two questions must be answered. One is 
that whether or not a thermal state can be estabhshed; another is that how fast 
O I thermalization of the system goes. The H-theorem ensures that the system evolves 
^ I into a thermal state as a result of 2-body elastic scatterings [1]. In a system where 
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3-body elastic scatterings are important, a transport equation needs a term that reflects 
the 3-body elastic scatterings. 
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where g is the color-spin degeneracy factor, g22 = n'^^v and g33 = riiJuout^- where n'^^^ 
(nont) is the number of identical final partons of 2-to-2 (3-to-3) scatterings and for 
the 3-to-3 scatterings is the number of identical initial partons except the parton in 
the distribution function fi. \ M.2-*2 P and | M.-^^^ P represent squared amplitudes of 
2-to-2 and 3-to-3 parton scatterings, respectively. 
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Define 

Hit) = J dfdpj{p,,f, t) In f{p,,r, t). (2) 

H{t) is negative for / < 1 and H{t) always decreases witli increasing time until it is 
independent of time because / obeys the transport equation. When H{t) does not rely 
on time, the distribution function takes a form of thermal distribution. Such reasoning 
proves that the system satisfying (1) is driven toward a thermal state by both 2-body 
and 3-body elastic scatterings. For a system that is governed by 2-body and 3-body 
elastic scatterings, the first question is answered. 

The answer to the second question depends on the squared amplitudes | J^2~>2 P 
and I A^3_^3 P which relate to dynamical processes. The very simple, instructive but 
unrealistic case is that the squared amplitudes are contants. Even in the case the 
integrals on the right-hand side of the transport equation are generally not zero. An 
initial distribution function is anisotropic in momentum space for a nonequilibrium 
system. If the two squared amplitudes are zero, the distribution function remains 
unchanged. Then thermalization is never possible. If the two squared amplitudes 
take small values, the distribution function changes slowly with increasing time and 
thermalization is slow. Large values of the two squared amplitudes will cause large 
change of the distribution function and the corresponding thermalization goes fast. 
Even though the squared amplitudes are constants, we still need to numerically solve 
the transport equation to get thermalization time which is the difference of the time 
when the initial distribution is given and the time when the equation produces an 
isotropic momentum distribution. Therefore, how fast a system thermalizes depends on 
dynamical processes like elastic scatterings and anisotropy of initial distribution. 

For a realistic case of interest, quark-gluon matter created in ultra-relativistic 
heavy-ion collisions, the squared amplitudes for 2-to-2 and 3-to-3 parton elastic 
scatterings have complicated dependences on parton momenta. Since the elastic 
scatterings take place among gluons, quarks and antiquarks with the same and/or 
different flavors, the study of how fast quark-gluon matter thermalizes becomes a very 
difficult task. However, from the discussion in the last paragraph, we immediately realize 
that thermalization of quark-gluon matter goes faster and faster if more and more types 
of parton elastic scatterings are included in | A42^2 P and | TVls^a p. A relevant change 
of thermalization time is shown by our recent works [21 [3l H] that have taken into 
account the four types of 3-to-3 elastic scatterings: quark-quark-quark, quark-quark- 
antiquark and quark-antiquark-antiquark for quark matter and gluon-gluon-gluon for 
gluon matter. The 3-to-3 elastic scatterings, thermalization time and summary are given 
in the next two sections. 

2. 3-to-3 Elastic scatterings 

The 3-to-3 quark scatterings at the lowest order are shown in figure 1 |2]. The left 
diagram indicates two-gluon-exchange induced scattering and the right involves triple- 
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gluon coupling. Initial (final) quarks may have the same or different flavors. While two 
or three initial (final) quarks have the same flavor, exchanges of quarks must be taken 
into account. 

Two typical processes of 3-to-3 quark-quark-antiquark scatterings at the tree level 
are shown in figure 2 [3]. A lot of processes are related to annihilation and creation of 
quark and antiquark. An exchange of the two initial (final) quarks needs to be considered 
if the quarks have the same flavor. Quark-antiquark-antiquark elastic scatterings are 
similar to quark-quark-antiquark elastic scatterings. 

Because of triple-gluon couplings and four-gluon couplings, the 3-to-3 gluon elastic 
scatterings are complicated. Four typical processes at the lowest order are exhibited in 
figure 3 [1]. Only the diagram B^^ involves two four-gluon couplings. The other three 
diagrams include one four-gluon coupling and two triple-gluon couplings. The diagram 
Bu has the four-gluon coupling related to two final gluons, B^+ to two final gluons and 
one initial gluon, and B 4(55) to one exchanged final gluon and two initial gluons. 

Except for the diagram B of which the squared amplitude is calculable by hand, 

we make fortran codes in Feynman gauge to derive the squared amplitudes of other 
diagrams that indicate the elastic scatterings of qqq qqq, qqq qqq, qqq qqq and 
ggg — ggg, respectively. While the squared amplitudes for the 3-to-3 elastic scatterings 
are used in the transport equation, the quark and gluon propagators are regularized 
with a screening mass [5] to remove Coulomb exchange divergence. 

3. Numerical solutions and summary 

Anisotropic parton momentum distributions are formed in initial central Au-Au 
collisions. Such anisotropy can be eliminated by elastic scatterings among partons. 
For a central Au-Au collision at ^/SNN = 200 GeV, we obtain from HIJING Monte 
Carlo simulation [6l [7] the initial gluon distribution that corresponds to a gluon number 
density of 19.4 fm"^. 



where Ra = 6.4 fm, go = 16 and the rapidity —5 < y < 5. Any of the up-quark, 
down-quark, up-antiquark and down- antiquark initial distribution functions is assumed 
to be one sixth of the gluon initial distribution function. Starting from the time tini, 
when anisotropic quark-gluon matter is formed and ending at the time tjso, when local 
momentum isotropy is established, the transport equation is solved. Local momentum 
isotropy, in other words, a thermal state is obtained if distribution function curves in 
different directions finally overlap through the thermalization time tiso — ^ini- For quark 
matter that is governed by the qq — > qq and qqq — > qqq elastic scatterings, thermalization 
time is about 1.8 fm/c; for quark matter interacted by antiquark matter via qq —>■ qq, 
qqq —>■ qqq and qqq — > qqq elastic scatterings, thermalization time is about 1.55 fm/c; 
for gluon matter the gg — >■ gg and ggg ggg elastic scatterings give a thermalization 
time of the order of 0.45 fm/c. While the qq — ^ qq and qqq qqq elastic scatterings give 
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a long thermalization time, the qq qq, qqq qqq and qqq — >• qqq elastic scatterings 
apparently shorten thermalization time of quark matter. 

In summary, the larger the squared amplitudes for the 2-to-2 and 3-to-3 elastic 
scatterings are, the faster thermalization goes. | A^a^s P for quark matter gets larger 
and larger if we do a successive inclusion of quark-quark-quark, quark- quark- ant iquark 
and quark-antiquark-antiquark elastic scatterings; correspondingly, thermalization time 
of quark matter becomes shorter and shorter. Rapid thermalization of gluon matter 
is obtained from the 2-to-2 and 3-to-3 gluon elastic scatterings. The study of 
thermalization leads to the importance of multi-parton elastic scatterings [8j in addition 
to the 2-body parton scatterings P [ini [HI [12] • We note here that we have not studied 
the contributions of qg ^ qg, qqg qqg, qqg — > qqg and qgg qgg since our strategy is 
that we first give an independent study of gluon matter thermalization or of quark matter 
thermalization and then study the influence of the interplay of quark matter, antiquark 
matter and gluon matter. Because the derivation of the squared amplitude for a 3-to-3 
elastic scattering is complicated and very time-consuming and numerically solving the 
transport equation needs months, we have to study the 3-to-3 elastic scatterings one 
by one. But the contributions of qg qg, qqg qqg, qqg — > qqg and qgg qgg 
are expected to shorten thermalization time of quark matter since they increase both 
I -^2^2 P and I Ais^s P, and the contributions will be considered in the coming years. 
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Figure 1. Quark-quark-quark elastic scatterings. 




Figure 2. Quark-quark-antiquark elastic scatterings. 
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Figure 3. Gluon-gluon-gluon elastic scatterings. 



